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This article describes the application of the digital signal processor 
as a voiceband signal processing element. The application chosen is 
one of the most stringent voice-frequency signal processing applica- 
tions in the telephone network— providing transmission treatment 
(e.g., gain, equalization, and echo control) for special service circuits. 
A detailed description of a prototype transmission treatment unit, 
which uses the dsp, is provided along with descriptions of the digital 
filter structures and filter synthesis techniques. Measured results for 
representative extreme cable facility cases are presented, showing 
that digital signal processing utilizing the dsp meets the telephone 
network transmission requirements for special service circuits. 

I. INTRODUCTION 

The first, and still most heavily used, transmission medium in the 
telecommunications network is copper wire. Twisted pair metallic 
cables of various gauges, lengths, and sizes make up the bulk of the 
loop plant and local exchange trunk plant. The application of electronic 
amplifiers or repeaters to provide gain to compensate for the atten- 
uation of signals on metallic cables was the first large-scale use of 
electronics in the telecommunications network. Repeaters are still 
used extensively. By necessity, the systems of the past that provided 
voice-frequency (vf) transmission treatment on a per-channel basis 
were analog systems. The large-scale use of sophisticated digital signal 
processing for these transmission treatment functions was precluded 
by cost, power consumption, and size of implementation. However, in 
recent years the capabilities of digital hardware have improved sig- 
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nificantly, primarily as a result of the rapid development of integrated 
circuit technology. 

The digital signal processor (dsp), 1 a vlsi device with a large number 
of logic circuits, is an example of the sophistication that is now possible 
with digital hardware. The inherent capabilities of the dsp have made 
it possible to consider the use of digital filters to replace some of the 
traditional analog network functions in vf transmission systems. This 
paper reports on the results of an experimental study in which the dsp 
is used for this purpose. Among the most desirable features of the dsp, 
and especially important for this study, is the ease with which it can 
be programmed under computer control to provide transmission treat- 
ment functions. 

The principal transmission signal processing functions are performed 
in vf repeaters for metallic transmission systems and in carrier ter- 
minal units (ctus) for transmission systems that contain a metallic-to- 
carrier interface. These systems provide signal processing for a variety 
of services. Typically, the most demanding signal processing require- 
ments derive from special services applications. 2 Special service circuits 
are engineered using all types of transmission media. Special services 
are a large and rapidly growing part of the telecommunications net- 
work. 

Performance objectives for a special service circuit are normally 
specified in terms of 1-kHz loss, attenuation distortion, echo distortion, 
and various other transmission parameters. If signal processing must 
be included in a circuit so that transmission objectives for the circuit 
can be achieved, the circuit is called a treated circuit. The prototype 
transmission treatment unit, described in this article, which uses the 
dsp to provide adjustable transmission treatment functions, is referred 
to as a digital treatment transmission unit (dttu). The description of 
the dttu and its capabilities begins in Section III. First, in Section II, 
an outline of treated circuits is presented since the performance goals 
for the dttu were based on the performance objectives for these 
circuits. 

II. INTRODUCTION TO VOICE-FREQUENCY CIRCUITS 

Table I lists some of the applications for metallic treatment systems. 
Figure 1 illustrates some of the possible circuit arrangements with 
transmission treatment. Each of the arrangements is assumed to be 
providing a special service circuit— in this case a foreign exchange 
trunk. Figure la shows a terminal repeater, with a switching system 
on one side and a 2-wire cable connected to a pbx on the other. Figure 
lb is smilar, but with the repeater placed at an intermediate location 
in the circuit. Figure lc shows the use of two repeaters in a long circuit. 

Finally, Fig. Id shows a circuit that contains both a carrier link and 
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Table I — Inserted Connection Loss (icl) 

objectives for typical 2-wire switched 

special services 



Switched Special Services 



ICL (dB) 



pbx-co trunk 
Foreign exchange trunk 
WATS trunk (to class 5 co) 
wats trunk (to class 4 co) 
Foreign exchange line 
wats line (to class 5 co) 
wats line (to class 4 co) 



3.5 
3.5 
3.5 
4.5 

3.5 

3.5 

4.5 



a 2-wire cable facility. This latter configuration, with a carrier link in 
tandem with a 2-wire metallic extension, is quite common. Carrier 
systems are designed to be essentially transparent for transmission 
purposes. Note that if the carrier link is removed in Fig. Id and the 
ctus retained, the resulting circuit arrangement is similar to that 
shown in Fig. la. 
The ctus that provide the 4- to 2-wire interfaces must provide 
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Fig. 1— Some example topologies for a foreign exchange trunk. 
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transmission treatment for the adjacent metallic cable. Carrier termi- 
nal units with analog treatment are currently available for both digital 
and analog carrier sytstems. A digital carrier system is an obvious 
candidate for digital transmission treatment since digital encoding and 
decoding of signals is already required. Digital switching systems, 
which have transmission characteristics similar to digital carriers, are 
also strong candidates for digital treatment of special service circuits. 
This article concentrates on the digital-to-analog 4- to 2-wire interface. 

2. 1 Transmission treatment objectives 

The success of the telecommunications network in providing satis- 
factory service places requirements or objectives on some of the basic 
electrical characteristics of the equipment used. The services listed in 
Table I must be engineered to have carefully controlled transmission 
properties. The main theme of this study is the use of the dsp to 
control the loss, attenuation distortion, and echo distortion of treated 
circuits. This is accomplished by implementing digital filters with the 
dsp, as is described in detail in the following sections. Although not 
considered here, the dsp could also be programmed to control delay 
distortion of a treated circuit. 

Of the services listed in Table I, the trunks have the most stringent 
performance objectives. A typical short haul treated special service 
trunk must be engineered to have a 1-kHz loss of 3.5 dB to an accuracy 
of approximately 0.5 dB. Furthermore, the loss of the circuit at 400 
and 2800 Hz, relative to the 1-kHz loss, should be within the following 
limits: at 400 Hz, -1.0 dB to 3.0 dB; at 2800 Hz, -1.0 dB to 4.5 dB. 
Additionally, the circuit loss should be relatively smooth between 
these two frequencies. The circuit should also be designed so that its 
loss response rolls off at the vf band edges to enhance stability margin 
and echo performance. For the configuration shown in Fig. la, crosstalk 
considerations restrict the gain of the repeater to a maximum of 6 dB 
at 1 kHz, and the loss of the cable to a maximum of 9 dB at 1 kHz. 

Since a 2-wire treatment unit is inherently a feedback device, some 
means must be provided to "balance" the unit, that is, significantly 
reduce the feedback, to ensure adequate stability margin and satisfac- 
tory echo performance for a treated circuit. In the dttu, a digital 
canceler network provides this balance function. Since gain added to 
a circuit amplifies an echo along with the desired signal, treated circuits 
must be better balanced than untreated circuits. Poor balance results 
in listener distortion on short circuits and talker echo on long circuits. 
At a 4- to 2-wire interface, the equipment should produce at least 15 
to 18 dB of loss between the 4-wire ports (e.g., ports 3 and 2 in Fig. 2), 
including the effects of gain, over most of the voiceband. The loss at 
the vf band edges may be somewhat less, although adequate stability 
margin must be maintained. 
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Fig. 2 — Balance objectives are achieved if there is at least 15- to 18-dB loss between 
ports 3 and 2 of the 4- to 2-wire transmission treatment unit. 



III. DIGITAL TREATMENT TRANSMISSION UNIT 

A block diagram of the dttu is shown in Fig. 3. The digital ports of 
the dttu could interface with a digital carrier or a digital switch, while 
the 2-wire analog port could interface with loop plant cable or trunk 
cable. The dttu consists of a dsp, a line interface unit (liu), and the 
appropriate logic devices to interface the dsp to the liu. All adjustable 
gain, equalization, and echo cancellation are provided by digital filters 
implemented with the dsp. The remaining dttu functions, such as 
encoding and decoding of signals, are provided in the liu. 

3. 1 Line interface unit transmission functions 

The experimental version of the liu contains a commercially avail- 
able 16-bit, full-linear digital-to-analog converter. Unpublished studies 
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Fig. 3— Block diagram of the digital treatment transmission unit (dttu). 
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by J. H. W. Unger have shown that the standard /x-255 encoding does 
not provide sufficient dynamic range for digital signal processing in 
some special service applications. The 16-bit linear encoding was 
chosen to avoid this problem. The converter is followed by a recon- 
struction filter that removes energy above 4 kHz. The reconstruction 
filter drives a fixed- gain analog amplifier that, in turn, drives the 2- 
wire cable through a 900-ohm transformer-coupled output. Addition- 
ally, the amplifier drives an analog compromise canceler, a device 
discussed in more detail in the next paragraph. Input signals from the 
cable are coupled through the transformer to a differential amplifier. 
The differential amplifier drives an antialiasing filter that removes 
energy above 4 kHz and also removes 60-Hz induction. The output of 
the antialiasing filter is sampled at 8 kHz and then converted to digital 
form by a 16-bit, full-linear analog-to-digital converter. 

The primary purpose of the compromise canceler is to provide some 
loss (approximately 6 dB or greater) between points "a" and "b" of 
the Liu for the universe of cables with which the dttu is expected to 
interface. Its transfer function is fixed, with one pole and no zeros. As 
compared to the dttu performance in the absence of the compromise 
canceler, the benefits obtained are twofold: 

(i) A large signal at point "a" is less likely to overload the analog- 
to-digital converter. 

{ii) The performance of the digital echo canceler in the dsp is 

enhanced. 

3.2 Digital signal processor transmission functions 

As mentioned above, all adjustable gain, equalization, and echo 
cancellation for the dttu are provided by digital filters in the dsp. 
These filters are: 

(i) an equalizer for the transmit direction of transmission, 
(ii) an equalizer for the receive direction of transmission, and 

(Hi) a canceler. 

The equalizers provide adjustable gain as well as adjustable equali- 
zation. The transfer functions of these three filters are dentoed by E,, 
E r , and C, respectively. These filters will be represented in many of 
the remaining figures in this article by the symbols shown in Fig. 4. 

IV. SELECTION OF FILTER FORMS FOR THE EQUALIZERS AND 
CANCELER 

The dttu must have the capability to provide transmission treat- 
ment for a large variety of cable facilities. Laboratory and computer 
simulations of a representative sample of treated transmission faculties 
have shown that the dttu has the required capability if the canceler 
is a 32-tap transversal filter and each equalizer is a fourth-order 
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Fig. 4 — Symbolic representation of dttu digital filters. 

recursive filter composed of two cascaded biquads. In Z-transform 
notation, the filter functions are: 



C= I c n z-", 

n-0 

P n FT am0r + a ^' Z ~ + a m2tZ~ 

*" ~ Ut " ~ i h — T 7 * h — ^2~' 

m=0 1 — O m \tZ — O m 2lZ 

JP /", TT flnl0r + a ™ lrZ _1 + °m2r2" 2 
Ar = G> J} — -rj — . 

m=0 A — O m \rZ — OmlrZ 



and 



(1) 

(2) 



Symbolic diagrams of these filters are shown in Figs. 5 and 6. (In the 
figures, x, y, Vi, and v 2 represent signal values, and are used for the 
discussion in Section V.) The equalizer pre-multipliers, G t and G r , 
allow flexibility in parceling out gain between the filter sections and in 




OUTPUT 

Fig. 5— Block diagram of a 32-tap transversal filter. 
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Fig. 6— Block diagram of two cascaded biquads. 

controlling coefficient magnitude, which is constrained by the dsp 
coefficient storage format to be less than 2.0. The coefficients in the 
three equations above are different for each different cable facility. 

A recursive structure, rather than a transversal structure, is used for 
the equalizers for one primary reason: in addition to providing gain 
and equalization in the frequency band 400 and 2800 Hz, for 2-wire 
transmission, it is also necessary that the equalizers provide significant 
attenuation at the vf band edges to enhance the stability margin 
performance of the dttu. This characteristic is more readily obtained 
with a low-order recursive structure, for the transmission system under 
consideration, than it is with a reasonable size transversal filter. 

A cascaded biquad structure, rather than a direct-form fourth-order 
structure, is used for the equalizers because the transfer function of 
the cascaded structure is less sensitive to coefficient quantization. In 
Refs. 3, 4, and 5 it is shown that if the unquantized pole locations or 
zero locations lie close to each other, coefficient quantization can cause 
a large shift in the pole or zero locations. If this shift occurs, the filter 
response will differ from the desired response. In the direct-form 
realization, the shift as a function of coefficient quantization is de- 
pendent on all the pole or zero locations. However, in the cascade 
form, the pole or zero shift in one section is independent of the pole or 
zero locations in the other sections. 

V. FILTER IMPLEMENTATION WITH THE DSP 

In this section, the implementation of the equalizers and canceler 
with the dsp is discussed. For a description of the DSP architecture see 
Ref. 1. There are four topics of interest with regard to the dsp program. 
These are: 
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(i) location of coefficients and delays in the ram, 
(ii) filter implementation, 
(Hi) coefficient loading, and 
(iv) program storage and execution time. 

5. 1 Digital signal processor RAM memory map 

Since the filter coefficients are different for each different cable 
facility, they must be stored in the dsp ram. Also, the delayed signal 
values must be stored in ram. Coefficients and delayed signal values 
are stored as shown in Fig. 7 so as to minimize the number of register 
sets required to access them. 

5.2 Filter implementation 

5.2.1 Equalizers 

In the dsp program used for the dttu, the three filtering operations 
are performed in the following order: 

(/) E r , 

(ii) C, 

(Hi) E t . 
The canceler is discussed below. In this section E r and E, are discussed. 
Since the program steps used to implement E r and E, are identical, a 
general description of the filtering steps in two cascaded biquads is 
presented. 

Figure 6 shows a block diagram of two cascaded biquads. The input 
signal is labeled x and the output signal is labeled y. Internal signal 
values are labeled v and u u with v (n - 1), v (n - 2), V\(n - 1), and 
Vi(n - 2) being delayed signals obtained from preceding filter opera- 
tions. The general sequence of operations used to filter x and obtain y 
requires nine additions and eleven multiplications for each value of x 
input to the filter. 

First, Uo(n) is obtained by forming the sum 

b 02 * v (n - 2) + boi * v (n - 1) + G * x(n) 

in the accumulator, rounding the result, and storing it in a temporary 
register. The rounded result is v (n). The products are performed in 
the product register before entering the accumulator. Then, vi(n) is 
obtained in an identical manner from the sum 

ao2 * v (n — 2) + aw * v (n - 1) + aoo * Vo(n) 

+ bi 2 * Vi(n - 2) + 6n * Ui(n - 1). 

Finally, y(n) is obtained from the rounded value of the sum 

ai2 * v\(n — 2) + an * Vi(n - 1) + aw * Vi(n). 
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♦SERVES SAME PURPOSE FOR CANCELER THAT G r AND G f SERVE 
FOR EQUALIZERS 

Fig. 7— Locations of the coefficients and delays in dsp ram. 



Before the next value of x enters the filter, v (n - 2) is updated by 
setting it equal to vo(n - 1); v (n - 1) is updated with v (n); Vi{n - 2) 
with vi{n - 1); and vAn - 1) with Vi{n). 

The output of the receive path equalizer is transmitted to the digital- 
to-analog converter and to the canceler input. 

5.2.2 Canceler 

Figure 5 shows a block diagram of a 32-tap transversal filter. The 
basic operation in a transversal filter is the multiplication of a tap 
coefficient by a delay output, followed by an addition in an accumu- 
lator. This operation can be accomplished with a sequence of instruc- 
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tions a = p + a p = axi * ryi;, where a and p are zero at the beginning 
of the sequence and rx and ry (coefficient and data address registers) 
are set to point to c 3i and x(n — 31), respectively. Because the dsp 
architecture does not permit a write to ram two instructions before 
reading a coefficient from ram, the delay shifts must be done separately 
from the filtering. 

Once the output of the canceler is computed, it is subtracted from 
the analog-to-digital output in ibuf. The contents of the accumulator 
are then transferred to the w register, ready for use by the transmit 
path equalizer. 

5.3 Coefficient loading 

A two-step procedure is used to load the dsp ram. In the first step, 
cO is set and the address of a ram location is transferred to the dsp. In 
the second step, cl is set and the coefficient is transferred and stored 
in the ram location that has the address transferred in step one. Each 
step takes one frame. The selection of the appropriate step is deter- 
mined by testing the cO and cl control bits. 

The procedure is repeated for all 54 coefficients. The coefficients are 
transmitted in order from address to address 53. Because the dsp 
ram is dynamic, and only one memory location is referenced every 
other frame, it must be refreshed every frame. This is accomplished by 
sequentially reading all locations. 

5.4 Program storage requirements and execution times 

The program to implement the filter functions and the coefficient 
loading routine requires 314 words of memory, with 21 memory loca- 
tions used for the receive path equalizer, 108 for the canceler, 36 for 
the transmit path equalizer, 85 for the coefficient loading and refresh, 
and 64 for miscellaneous operations (pc sets, no-ops, etc.). In the 
experimental version of the dttu, code was written to be clear rather 
than efficient in memory requirement. There are some techniques that 
in a final product could reduce the amount of memory required. For 
example, in doing the canceler tap update, the instruction rdi = ryp; 
was used. This requires two words of storage per "auxiliary" instruc- 
tion. By using rdi = ryp a = p p = olx * c; the amount of storage for 
the update sequence is cut in half, since this is a "normal" instruction. 
The a = p p = olx * c; is just a "fill" and accomplishes no useful 
operation. 

The filtering portion of the program requires 106.40 /is, with 12.8 /is 
for the receive path equalizer, 56 /is for the canceler, 12.8 /is for the 
transmit path equalizer, and 24.8 /is for miscellaneous operations. This 
leaves 18.6 /is for other use. One application might be a self-diagnostic 
for the chip (e.g., check ram or the au). 
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VI. LOSS DEFINITIONS 

A special service circuit must provide a high-quality channel for vf 
signal transmission. Loss measurements are normally used to assess 
the performance characteristics of such facilities since many of the 
performance requirements (or objectives) for the facilities are specified 
in terms of loss, e.g., return loss (or balance and echo performance), 1- 
kHz loss, and attenuation distortion. For an all-analog system, well- 
known definitions and measurement procedures exist for each of these 
loss types. A service provisioned with the dttu is not an all-analog 
system, but is, instead, a mixed analog-digital system. Since measured 
results of the characteristics of this system are important for all of the 
remaining material in this article, some of the conventions and defi- 
nitions used for loss measurements in a mixed analog-digital system 
are presented in this section to aid the reader in interpreting the 
measured results. 

Two additional, and related, subjects are also discussed in this 
section. These are: Measurements required to determine filter coeffi- 
cients and dttu loss scaling. Loss scaling (signal level control) is 
necessary for the dttu since a digital processing system has a limited 
dynamic range (determined by the 16-bit accuracy of the analog-to- 
digital and digital-to-analog converters, in this case). 

A simplified topology of a typical 2-wire treated facility, with treat- 
ment provided by the dttu, is shown in Fig. 8. The dttu ports are 
labeled 1 (2-wire port), 2 (transmit port), and 3 (receive port), with 2 
and 3 being digital ports and 1 an analog port. Depending on the 
service being provided, the terminating equipment at the far end of 
the cable (port 4) could be a central office switch or various types of 
customer premises equipment (telset, pbx, etc.). 

6.1 Loss in a mixed analog-digital system 

The results of three loss measurements, denoted by L,, L r , and L c , 
are necessary for characterizing the transmission performance of the 




Fig. 8— Simplified topology of a typical 2-wire treated facility. 
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Fig. 9 — Transmission path and facility configuration used for measurement of L,. 

facility shown in Fig. 8. The transmission paths associated with L t , L r , 
and L c are shown in Figs. 9, 10, and 11, where, as the figures show, 

L t = loss from far end of cable to transmit port of the 
dttu (analog-digital loss), 

L r = loss from receive port of the dttu to far end of 

cable (digital-analog loss), and 

L c = loss from receive port of the dttu to the transmit port 
of the dttu (digital-digital loss). 

To define L,, let a 0-dBm analog sine- wave generator of internal 
resistance R be attached to the far end of the cable. Internal resistance 
R is either 900 or 600 ohms, depending on the impedance characteristic 
of the terminating equipment. At port 2, a sampled representation of 
this sine-wave will appear, with amplitude and phase characteristics 
determined by the transfer function of the cable-DTTU system. The 
information required for determining L t is the peak of the sine-wave 
signal at port 2. It is unlikely that the peak of this signal can or will be 




Fig. 10 — Transmission path and facility configuration used for measurement of L r . 
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Fig. 11— Transmission path and facility configuration used for measurement of L,, 

encoded; therefore, the samples, which are considered to have integer 
representations, must be analytically processed to determine the peak. 
If the real number representation of the signal peak at port 2, as 
analytically determined, is called V, then L, is defined as 

L t = 20 log(VV V) - 20 log(V>./Vo) = 20 \o%{V Q /V), 

where V is the level representative of dBm at port 2 and V fs is the 
full-scale level (32767.0 for a 16-bit system); V is less than V fa , and for 
the system under consideration is defined by the relationship 

201og(V /s /Vo) = 3dB. 

This definition applies to both the transmit and receive ports and is 
consistent with the convention used in digital carrier systems where 
the code is jn-255. 

To define L r , let a digital sine-wave of peak level V be input at port 
3. Then, an analog sine-wave will appear at the end of the cable. Its 
power in dBm, measured across resistance R, is L r . 

Finally, to define L c , let a digital sine-wave of peak level Vz be input 
at port 3. Then, a digital sine-wave of the same frequency will appear 
at port 2. If the peak level of this signal is V2, then 

L c = 201og(V 3 /V 2 ). 

When a measurement of L c is made to assess the echo and stability 
performance of the transmission facility, the cable will normally be 
terminated in a standard termination that is representative of the 
impedance of the terminating equipment at the far end of the cable. 

6.2 Measurements required to determine filter coefficients 

To set the equalizer and canceler filters for treatment of a metallic 
cable pair, the transmission characteristics (loss and input impedance 
vs. frequency) of the facility must be known. These may be determined 
either from a mathematical model of the cable or from measurements 
of its characteristics. A typical mode of operation for a facility to be 
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provisioned with a dttu would be one in which the characteristics of 
the cable are measured. Since the Liu is interposed between the dsp 
and the cable facility, its transmission characteristics must be included 
in the measurements. 

In the laboratory arrangement used for this study (and in an envi- 
sioned practical application), it is most convenient to make these 
measurements, which are identical to the measurements discussed in 
Section VI, utilizing digital access points in the dsp. Since measure- 
ment access is through the dsp, the equalizers are set to unity and the 
canceler is set to zero. For this condition, the three loss measurement 
results are denoted by Lto, Lro, and Leo- 

6.3 Digital treatment transmission loss scaling 

A transmission treatment unit is required to operate as a low-noise, 
linear amplifier for a wide range of signals expected to be flowing in 
the telephone network. Therefore, signal levels in the unit must be 
controlled so that a large signal is not over-amplified and hence 
distorted and a small signal is not overly attenuated toward the noise 
floor of the unit before being amplified, thereby degrading its s/n. In 
the dsp, one of the most critical interfaces where signal levels must be 
controlled is the dsp-liu interface. A large signal incident on the 
receive port should not be over-amplified by E r and therefore overflow 
and be distorted at the dsp-liu interface. Additionally, a small signal 
incident on the 2-wire line port should not be overly attenuated in the 
Liu before being encoded and subsequently amplified by E,. 

To control the signal levels at the dsp-liu interface, the dttu is 
loss-scaled, i.e., signal loss through the unit is adjusted to achieve the 
desired level control. Fixed-loss (or gain) amplifiers in the liu are used 
for this purpose. 

The magnitude of the loss scaling used for the experimental version 
of the dttu can be exhibited by demonstrating the effect it has on loss 
measurements performed on a cable-DTTU system when the equalizers 
are set to unity. In the frequency range 200 to 3400 Hz, where the 
antialiasing and reconstruction filters introduce negligible frequency 
shaping, 

Lro = cable loss — 7 dB 

and 

Lto = cable loss + 3 dB. 

Between 3400 and 4000 Hz, I& and Lto remain approximately 10 dB 
apart; however, filter attenuation contributes significantly to their 
specific values. 
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VII. EQUALIZER AND CANCELER COEFFICIENTS 

This section discusses the procedures for determining equalizer and 
canceler coefficients from the measured facility transmission charac- 
teristics. These procedures are based on minimum-least-square curve 
fit techniques. 

7. 1 Overview of the curve fit procedures 

The canceler and equalizer coefficients are determined by separate, 
frequency-domain, minimization algorithms. The goal, in each case, is 
to minimize a penalty function. These penalty functions are: 

Pc= I \C(F) - T C (F)\ 2 , (3) 

P,= I (\E l (F)\-\T l (F)\) 2 , (4) 

in* 

Pr= I (\Er(F)\-\Tr(F)\) 2 , (5) 

(F) e 

where the target functions T c , T,, and TV are dependent on the 
transmission characteristics of the cable-DTTU system. Functions T, 
and T r are also dependent on the attenuation distortion objectives and 
1-kHz loss objective for the facility that is to receive treatment. {F} c 
denotes a set of frequencies for which P c is to be minimized, and [F) e 
denotes a set of frequencies for which P t and P r are to be minimized. 
Extensive study of the results achieved from minimization of the 
penalty functions listed above, for various frequency sets and for 
various cables, has shown that satisfactory curve fits can be obtained 
if 

{F)c = 100 to 3900 Hz in 100-Hz increments 
and 

{F} e = 100 to 3700 Hz in 300-Hz increments. 

In eqs. 3, 4, and 5, the canceler and equalizer functions are expressed 
in the frequency domain, i.e., the transform variable z has been 
restricted to the unit circle and is equal to 

exp(J2^F/F s ), 

where F s is the sampling frequency (8 kHz). 

As is shown below, minimization of P c results in a linear solution for 
the coefficients of C, while equalizer coefficients must be determined 
with a gradient search technique. Additionally, it is shown that only P, 
or P r need be minimized since, for a 2-wire facility, E t and E r differ 
only by a constant gain factor which is a result of the loss scaling used 
for the dttu. 
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7.2 Canceler target function 

In both magnitude and phase, T c is equal to the transfer function 
from the receive port to the transmit port of the cable-DTTU system 
when the equalizers are set to unity and the canceler is set to zero. 
Using the notation introduced in Section VI, the magnitude of T c is 

\T C \ = 10~ Lco/20 . 

Figure 12 shows plots of T c , in dB's, for two example cable cases and 
in Figs. 13a and 13b the corresponding phase plots are shown. 

7.3 Equalizer target functions 

Functions T t and T r are determined from the following information 
for the facility that is to receiver treatment: 
(i) 1-kHz loss objective, 
(ii) attenuation distortion (ad objectives for the frequency set 

{F}e, 



2.16/iF 




-30 



CABLE 1 = 18 kft OF 26-GAUGE NONLOADED CABLE 
CABLE 2 = 36 kft OF 24H88 LOADED CABLE 



Fig. 12 — T, versus frequency for two cable cases. 
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Fig. 13a— Phase plot of T c for cable 1. 
Fig. 13b— Phase plot of T c for cable 2. 



(Hi) L/o and L^ (see Section VI) for the frequencies from 400 to 2800 
Hz in 300-Hz increments. 
However, knowing T t is equivalent to knowing T r since 

(i) for all frequencies in the voiceband, for a treated 2-wire facility, 
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the transmit direction loss of the facility should equal the receive 
direction loss (measured on an end-to-end basis), 

(ii) Lto and Lro have the same shapes above 200 and below 4000 Hz, 
but differ in absolute magnitude by a value that is independent of 
frequency or cable type (see Section VI). 

Therefore, T, and TV have the same shapes but differ by a constant. 
The same conclusion can therefore be drawn about E t and E r . 

Since E r can be determined from knowledge of E, and vice versa, 
the notation distinction between E t and E r is dropped for much of the 
remaining discussion. The following simplified notation will be used: 
(i) E = equalizer function, 

(ii) L = loss data used in determining equalizer coefficients (i.e., L«> 
or Lm), 

(Hi) T = equalizer target function, and 

(iv) P = equalizer penalty function. 
Therefore, 



P= I (\E(F)\-\T(F)\) 2 , 



where 



e = g n 



I _1 o 

a m a + a, m \Z + O-miZ 



m-0 1 — Oml Z — UmzZ 

And finally, T will most often be expressed in dB's with the notation 

D= 20 log | T\. 

Now the discussion returns to the main topic, the determination of T 
(i.e., D). First, the facility objectives are considered. 

For 2-wire vf treated services, the loss objective is specified at 1 
kHz and the ad objectives are specified at 400 and 2800 Hz. The ad 
objectives offer a "window" of allowed values. Since the procedure for 
determining equalizer coefficients is a curve fit, specific ad objectives 
must be chosen. The choices are: 

1 - dB roll-off at 400 Hz 

2 - dB roll-off at 2800 Hz. 

Additionally, to ensure that the equalizers are well-behaved between 
400 and 2800 Hz, it is necessary to specify ad objectives at other 
frequencies between these endpoint frequencies. The result is an ad 
objectives "curve" specified at 300 Hz increments between 400 and 
2800 Hz. The ad objective at 1 kHz is dB. The ad objective curve 
used for this study, filled-in between the 300-Hz increments, is shown 
in Fig. 14. For the remainder of the discussion, the curve depicted in 
Fig. 14 is called 0„. 
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Fig. 14 — Attenuation distortion objective curve 0„. 



The next quantity needed for determining D is L, which is to be 
obtained by measurement. It is required that L be known for the set 
of frequencies from 400 to 2800 Hz in 300-Hz increments. Figure 15 
shows plots of L, normalized to dB at 1 kHz, for two different cable 
facilities. 







CABLE 1 = 18 kft OF 26-GAUGE NONLOADED CABLE 
CABLE 2 = 36 kfl OF 24H88 LOADED CABLE 



0.4 



0.8 1.2 1.6 2.0 

FREQUENCY IN KILOHERTZ 



Fig. 15 — L versus frequency for two example cables. 
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The quantities discussed up to this point are sufficient for determin- 
ing D in the frequency band from 400 to 2800 Hz. In this frequency 
band, 

D = L-O a -K h 

where K x is the 1-kHz loss objective for the treated circuit under 
consideration. 

Outside the frequency band from 400 to 2800 Hz, a shape for D is 
chosen that forces the equalizers to roll off at the vf band edges. 
Figure 16 shows D curves, normalized to dB at 1 kHz, for two 
different cable facilities. 

7.4 Algorithm for minimizing P c 
The goal of the minimization algorithm for P ( is to obtain a set of 
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Fig. 16 — D = 20 Log | T\ versus frequency for two example cables. 
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canceler coefficients that will result in maximum loss from port 3 to 
port 2 of the dttu. The transfer function from port 3 to port 2 is 

-i=£ 2 (T,-C). 

V3 

Therefore, the function that must be minimized, in a least squares 
sense, is 

^•= I \E\ 2 \T C -C\\ 

It is desirable that C not depend on the transfer functions of the 
equalizers; therefore, all equalizers are set to unity. The resulting 
function that must be minimized is 



Pc- 


« Z \T C -C\ 2 

lF) e 


To minimize P c , set 




^ = 0, 

dC m 


m = 0, 1, 2, • 


The result is a matrix equation 




BC = D, 


where the elements of B are 





31. 



b mn = £ cos|>F(m - n)/F s ], 

C is the vector of tap weights, and D is a vector with elements 
d m = I ne[T c (F)exp(j27rFm/F s )]. 

The solution for C is 

C^B-'D. 

7.5 Algorithm for minimizing P 

For vF-treated services, phase equalization is usually not required; 
therefore, it is necessary only to obtain a best fit for the magnitude of 
E. In a least squares sense, the function that must be minimized is 

P= I (|£(F,)|-|r(F,)|) 2 , 

where F, = (300i — 200) Hz. The simple method used for minimizing 
P c does not yield a linear system of equations when applied to P. 
Instead, minimization of P requires a nonlinear optimization proce- 
dure. The procedure chosen was developed by Steiglitz. 6 Steiglitz's 
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procedure uses the Fletcher-Powell' optimization technique and a 
novel approach for choosing initial conditions. 

To apply this procedure, the equalizer transfer function is first 
rewritten in the form 

iP/ir a v, a * + *° 2 ~' + XlZ ~ 2 * + XaZ ~ X + X5Z ~ 2 AWtW 1t\ ia\ 
E(F,A,X)=A- n =9- =1 — 2 =Atl(b,X), (6) 

1 - XtZ - Xaz z 1 - XbZ x - XiZ * 

where 

X = ( x , x h x 2 , x 3 , x 4 , x 5 , x 6 , x 7 y. 

(The prime denotes transpose.) P can now be written in the form 

P= £ (\AH(F i ,X)\-\T(F i )\) 2 . 
1=1 

Since |i4| appears in a linear fashion, it is easy to show that the 
optimum value of | A \ is 

13 



I \H(F„X)\\T(F,)\ 



IAI=- 



Z \H(F,X) 



!=i 



The sign of A is irrelevant and is therefore chosen to be positive. Since 
A can be precisely determined from knowledge of X and the 1-kHz 
gain of the equalizer, the nonlinear optimization algorithm is used only 
to determine X. 

For most nonlinear optimization techniques, the primary difficulty 
is in choosing an initial set of coefficients that guarantees convergence 
in a reasonable number of iterations. The following method for choos- 
ing initial coefficients has worked quite well for minimization of P. 

The first step is to find an optimum solution for the simpler equalizer 
function 

— 1 — 2 

El{Ft Alt Y)=A i ] +y ° Z . 1 +yiZ . 2 = AtfAF, Y), 
1 - yiz — y$z 

where 

Y = (jy , yx, J2, ys)'. 

The initial values of the coefficients, Y l , are chosen to be 

Y 1 = (0.0, 0.0, 0.0, 0.0)'. 

The final values of the coefficients, obtained by minimizing 

Pi(A h Y)=l (\AiHi(F h Y)\-\T(F i )\) 2 



1=1 



with the same algorithm that will be used to minimize P, are denoted 

by 
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Y* = (yty*iytyty. 

Following the minimization of Pi, P is minimized, with the initial 
coefficients being 

X 1 = (yt, yt, yt, yt, 0.0, 0.0, 0.0, 0.0)'. 

The final value of X is denoted X*. 

Figure 17 shows a flow chart of the optimization algorithm used to 




ENTER FUNCTION FP. 
CURRENT VALUES OF 
PANDX RETURNED. 



ENTER FUNCTION PZ 






|n 


Pq 


= p 


1 


ITER = 


ITER+ 1 




P = PENALTY 
FUNCTION 
FOR PRESENT 
INTERATION 
P = PENALTY 
FUNCTION 
FOR PREVIOUS 
INTERATION 
„ N = MINIMUM 
PENALTY 
FUNCTION 
FOR ALMOST 
PERFECT 
SOLUTION 



Fig. 17 — Flowchart of algorithm used to minimize P, and P. 
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minimize Pi and P. ncoefs denotes the number of coefficients that 
are being determined for any one pass through the algorithm, with 
ncoefs = 4 for minimization of Pi and ncoefs = 8 for minimization 
of P. In the algorithm, the common notation for both penalty functions 
is P. The reader should note the calls to functions FP and PZ in the 
algorithm. 

Function FP contains the Fletcher- Powell procedure. The FP re- 
turns the current value of P and the corresponding array of coefficients. 
Within FP there are a maximum of 25 updates to the coefficients. The 
array X is used to store the array Y when ncoefs = 4. 

The purpose of function PZ is to invert all poles and zeros that are 
outside the unit circle to the inside of the unit circle and to move all 
poles and zeros that are on the unit circle to a small distance inside. If 
any poles or zeros are outside the unit circle, the actions taken in PZ 
changes the phase of E, but not the magnitude, which is the important 
quantity in this case. 

The primary justification for the algorithm just described is that it 
has never failed to yield a suitable minimum for P, for realistic target 
functions. The algorithm takes 10 to 30 seconds of CPU time, depending 
on the cable facility, on a DEC PDP 11/70, and is written in the C 
language. 8 

The output of the algorithm, X*, does not contain enough informa- 
tion for a direct determination of the optimum set of equalizer coeffi- 
cients. The information contained in the vector X* is sufficient only 
for a determination of the poles and zeros of E. Additionally, the 
algorithm has not necessarily yielded an ordering of the poles and 
zeros that is best suited for achieving the noise and distortion objec- 
tives for the facility. These issues are now considered. The outcome is 
equalizer transfer functions in their final form, i.e., the form expected 
by the dsp program. 

7.6 Procedure for putting the equalizer functions in their final form 

When the 1-kHz gain of the equalizers is determined, E t and E r can 
be expressed in the form 

Et = A,H(X*), (7) 

E r - A r H(X*), (8) 

where H was defined in eq. 6. The object of the discussion in this 
section is to outline a procedure for putting eqs. 7 and 8 in the form of 
eqs. 1 and 2. For the types of signals expected to be flowing through 
the equalizers, the goals of this procedure are as follows: 

(i) There should be minimum degradation of the s/n of the signals. 
(ii) There should be a minimum amount of signal clipping at the 
critical internal nodes (to be defined below) of the equalizers. 
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These two goals can operate at odds with each other since the 
second goal can often be achieved at the expense of the first. In 
arriving at the filter forms given in eqs. 1 and 2, the magnitude of each 
coefficient must be constrained to be less than 2.0. 

The procedure discussed below does not necessarily have general 
applicability; instead, it is intended for the specific application de- 
scribed in this article. The justification for the procedure is that it 
yields the intended goals, as has been born out by experience. 

The first stage in the procedure is to select an ordering of the poles 
and zeros of the equalizer. It should be noted that eqs. 7 and 8 can be 
expressed in the form 

II (1 - z n z~ l ) [1 Ni 

Et = A,^- = A, "'" 



4 —' 1 

n (i - Pnz~ i ) n & 

71=1 J=0 



i-0 

where No, Ni, Do, D\ are second order functions of 2. If all the z n 's are 
real, there are six ways of forming Ni, while if at least two of the z„'s 
are complex, there are only two ways of forming Ni. Obviously, the 
same statements apply to Z), as to the p n 's. 

To select the ordering of the poles and zeros of E, and E r , consider- 
ation must be given to the transfer functions from the input of the 
filter to the critical internal nodes 9 of the filter. For two cascaded 
biquads, there are two critical internal nodes, located at the output of 
the accumulator and preceding a multiplier (see Fig. 18). At the critical 
internal nodes, overflow must be prevented for large signals input to 
the filter, and over-attenuation must be prevented for small signals 
input to the filter. The reader should note that these goals are identical 
to those for dttu loss scaling described in Section VI. The transfer 
functions to the critical internal nodes of two cascaded biquads are 
related to the N's and D's as follows: 

Ti is proportional to \/D\, 

Ti is proportional to N\/(D\D'i). 

For the system under consideration, these transfer functions are picked 
according to the following criterion. Each critical node transfer func- 
tion is chosen such that, in the frequency range 200 to 3400 Hz, the 
difference between the largest value of the magnitude of the transfer 
function (in dB's) and the smallest value of the transfer function is as 
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Fig. 18 — Critical internal nodes of two cascaded biquads. 

small as possible for the choices that are available for the ordering of 
poles and zeros (i.e., the response is as flat as possible under the 
constraints imposed). 

For either T\ or T2 there may be six choices or two choices, 
depending on the number of complex poles and zeros. The rationale 
for the criterion given above is that it minimizes the possibility that 
gain- scaling of the filter sections will result in overflow or a degradation 
in s/n at any of the frequencies in the voiceband. 

When the ordering of poles and zeros has been completed, the next 
step is to parcel out the gain factors A t and A r . The notation used for 
the equalizer functions is 



a m ot + a m uz ' + a m 2tz 



E r =G r 



Do D* 

GorNo GlrN, 

Do Di 



1 — bmitZ — b m 2fZ 

-1 1 „ --2 



m=0 

_L a m or + a m \ r z~' + a m 2rZ 

= G r 11 



m=0 \ — b m \ r Z — b m 2rZ 



where 



GtGotGu = At, 

GrGorGlr = A r . 

For the two directions of transmission, different criteria are used to 
gain-scale the filter sections. 
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The signals incident on E, are small, having been attenuated by the 
cable. Therefore, to ensure that s/n for these signals is not degraded 
further, G t and Got are chosen at their maximum allowed values, 
determined by the coefficient magnitude constaint. Then, having 
picked Gt and G 0I , Gu is 

Ax 



G„- 



GtGo 



The signals incident on E r are large compared to those incident on E,; 
therefore, to ensure a minimum amount of clipping at the critical 
internal nodes of the filter, G r and Go r are chosen to be as small as 
possible and G\ r is chosen at its maximum allowed value. G ir is chosen 
first. Then G r and Gor are chosen as follows: 

(9) 

unless, of course, the resulting Go r is too large (which is very unlikely) 
to allow GorNo to satisfy the coefficient magnitude constraint. If G 0r is 
too large, as defined by eq. 9, it is chosen at its maximum value, with 
G r chosen as 




GorGi 



At this point enough information is available for determination of the 
final values of the equalizer coefficients as well as the pre-multipliers 
G t and G r . The next important topic is the transmission performance 
of facilities for which treatment is provided by the dttu. 

VIII. TRANSMISSION TREATMENT CAPABILITIES OF THE DTTU 

In this section the capabilities of the dttu in providing transmission 
treatment for metallic facilities is presented. Four cable cases are used 
in the examples. These cables were chosen because they have charac- 
teristics, either in loss or impedance, that are representative of "worst 
case" cables allowed in treated services. The first two cable cases were 
used for the example measurement results for the discussion in Section 
VII. For all cases, the cables were simulated in the laboratory with 
artificial cable kits. 

8.1 Equalization 

Figs. 19 through 22 show the results of loss measurements performed 
on a cable-DTTU system for each of the four cables. Each figure 
contains two plots, labeled "a" and "b". These plots are of 
(i) unequalized receive direction loss (a) and 
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(ii) equalized receive direction loss (b). 
For each cable, the equalized loss for the transmit direction is identical 
to the equalized loss for the receive direction, except below 200 Hz 
where the transmit direction loss is greater because of the low-fre- 
quency roll-off of the anti-aliasing filter. The unequalized transmit 
direction loss is approximately 10 dB greater than curve "a" for 
frequencies above 200 Hz. 

For each curve labeled "6", the roll-off at 400 and 2800 Hz is 
displayed. As is evident, the roll-off results are close to the 1-dB 
objective at 400 Hz and the 2-dB objective at 2800 Hz (see Section 
VII). The 1-kHz loss objective for each facility is 3.5 dB. In the 
laboratory, the high-quality equalized performance shown for the four 
example cables has been consistently achieved for a large selection of 
other cables. 



8.2 Echo cancelation and balance 

Figures 23 through 26 each show the results of three separate 
measurements of loss from port 3 to port 2 of the dttu, with each of 
the example cables in turn attached to the 2-wire port of the unit. For 
all measurements, the cable termination is 600 ohms in series with 2.16 
uF, one of the standard terminations used for vf transmission systems. 
For the measurement results labeled "a", the equalizers were set to 
unity and the canceler was set to zero, leaving only the effect of the 
fixed compromise canceler. For the measurement results labeled "6", 
the equalizers were also set to unity, but the canceler was set to the 
result obtained through the minimization algorithm discussed in Sec- 
tion VII. For the measurement results labeled "c", equalizers and 
canceler were set to the results obtained through the minimization 
algorithms. 

In Figs. 23 through 26, note that, on the whole, the effectiveness of 
the canceler increases with frequency (see curves "6"). This type of 
response counteracts the gain of the equalizers which also increases 
with frequency for the nonloaded cables. The low-frequency effective- 
ness of the canceler could be improved by increasing the number of 
taps. However, increasing the number of taps is not necessary since 
the roll-off in the equalizers at these low frequencies yields the desired 
improvement in echo performance. The curves labeled "c" illustrate 
the echo performance of the facility when full treatment is applied. 
Obviously, the 15- to 18-dB echo objective (see Section II) is easily 
met for all facilities. The excellent performance shown for the example 
cables has been achieved for a large number of cables studied in the 
laboratory. 
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Fig. 23 — Loss measurements from receive port to transmit port (cable 1). Plot (a) 
represents equalizers set to unity, canceler to zero. Plot (b) represents equalizers set to 
unity, canceler set to result obtained with miriimization algorithm. Plot (c) represents 
equalizers and canceler set to results obtained with minimization algorithms. 
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Fig. 24 — Loss measurements from receive port to transmit port (cable 2). Plot (a) 
represents equalizers set to unity, canceler to zero. Plot (b) represents equalizers set to 
unity, canceler set to result obtained with minimization algorithm. Plot (c) represents 
equalizers and canceler set to results obtained with minimization algorithms. 
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Fig. 25 — Loss measurements from receive port to transmit port (cable 3). Plot (a) 
represents equalizers set to unity, canceler to zero. Plot (b) represents equalizers set to 
unity, canceler set to result obtained with minimization algorithm. Plot (c) represents 
equalizers and canceler set to results obtained with minimization algorithms. 



IX. CONCLUSIONS 

It is clear from the studies reported in this article that the dsp is 
capable of providing gain, equalization, and echo canceling for vf 
special services. Techniques have been developed for digital filter 
synthesis and implementation and for incorporating the dsp into a 
practical system. Applications to specific systems depend only on the 
issues of cost and power consumption relative to other techniques. 

X. ACKNOWLEDGMENTS 

The authors wish to thank the following people for the valuable 
contributions they made to this project: D. C. Watkins, R. J. Sanfer- 
rare, R. J. Gallant, R. L. Overstreet, J. Aagesen, M. R. Aaron, J. R. 
Boddie, and S. M. Walters. Special thanks go to D. C. Watkins and R. 



VOICE-FREQUENCY TRANSMISSION 1617 




2.16 /(F 



0.5 



1.0 



1.5 2.0 2.5 

FREQUENCY IN KILOHERTZ 



Fig. 26 — Loss measurements from receive port to transmit port (cable 4). Plot (a) 
represents equalizers set to unity, canceler to zero. Plot (b) represents equalizers set to 
unity, canceler set to result obtained with minimization algorithm. Plot (c) represents 
equalizers and canceler set to results obtained with minimization algorithms. 
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